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Abstract: The Pacific abalone, Haliotis discus hannai, is a popular food in Eastern Asia. 
Aquacultural production of this species has increased because of recent resource declines, 
the growing consumption, and ongoing government-operated stock release programs. 
Therefore, the genetic characterization of hatchery populations is necessary to maintain the 
genetic diversity of this species and to develop more effective aquaculture practices. We 
analyzed the genetic structures of five cultured populations in Korea using six microsatellite 
markers. The number of alleles per locus ranged from 15 to 64, with an average of 23.5. The 
mean observed and expected heterozygosities were 0.797 and 0.904, respectively. The 
inbreeding coefficient Fis ranged from 0.054 to 0.184 (mean Fis = 0.121 ± 0.056). The 
genetic differentiation across all populations was low but significant (overall F ST = 0.009, 
P < 0.01). Pairwise multilocus Fst tests, estimates of genetic distance, and phylogenetic and 
principal component analyses did not show a consistent relationship between geographic and 
genetic distances. These results could reflect extensive aquaculture, the exchange of breeds 
and eggs between hatcheries and/or genetic drift due to intensive breeding practices. Thus, 
for optimal resource management, the genetic variation of hatchery stocks should be 
monitored and inbreeding controlled within the abalone stocks that are being released every 
year. This genetic information will be useful for the management of both H. discus hannai 
fisheries and the aquaculture industry. 
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1. Introduction 

The abalone is a marine gastropod of the genus Haliotis. Nearly 100 Haliotis species inhabit both 
tropical and temperate waters worldwide [1]. Among the more than 30 species of commercially 
available abalone, the Pacific abalone {Haliotis discus hannai), which is distributed in the coastal areas 
of Eastern Asia, including Korea, Japan, and China, is one of the most popular species [2] . In response to 
a decrease in natural resources in the face of increasing demand, the aquacultural production of abalones 
has increased [3]. Culturing of H. discus hannai on a mass scale in Korea began in the 1980s, and has 
rapidly developed in the past decade. The majority of the annual yield now comes from aquaculture. The 
annual production of cultured H. discus hannai reached approximately 7500 tons in 2009 [4]. In addition 
to increased abalone farming, the Korean government has sponsored artificial abalone seed release for 
coastal abalone resource enhancement since 1997. More than five million Pacific abalone seeds reared in 
hatcheries were released into the Korean coastal sea areas in 2007 [5]. 

Natural populations of H. discus hannai have been overexploited to the extent that wild Pacific 
abalone is difficult to obtain in Korea. Pacific abalone seed is usually produced in hatcheries using 
reared adults as broodstock. Hence, hatchery production of abalone raises concerns regarding the 
maintenance of genetic diversity among cultured stocks, especially because their seedlings are released 
into natural habitats and thus could potentially alter the genetic structure of natural populations [6]. 
Despite a long history of aquaculture in Korea, the genetic diversity of hatchery stocks remains 
unknown. Therefore, an investigation of genetic variation in cultured abalone stocks is urgently needed 
for successful hatchery management, the production of high-quality abalone and to avoid reductions to 
the genetic variation present in aquaculture stocks. 

The monitoring of genetic variation among marine resources, especially in species for which artificial 
stocks produced by aquaculture are used for the restoration of natural resources, is essential to ensure 
that stock enhancement programs successfully preserve genetic diversity. This monitoring necessitates 
the development of genetic markers that can be used to assess genetic variation among populations and 
to prevent the release of mixed hatchery seeds, which cannot be detected by eye. 

Among available genetic markers, microsatellites are recognized as an essential tool in population 
studies because of their useful properties, such as high levels of polymorphism, codominant inheritance, 
and good reproducibility [7,8]. Over the past decade, microsatellites have produced promising results in 
studies of genetic variation in many marine species [9-13]. 

Until recently, microsatellite markers have been developed in the Pacific abalone [14-18], and 
genetic variability of hatchery stocks in Pacific abalone has been analyzed [6,17-20]. Microsatellite 
markers could sensitively detect the reductions of genetic variability on allelic diversity and mean 
heterozygosity. Highly significant Fst values were observed between hatchery stocks in Japan and 
China [6,19,20]. There are only two examinations for genetic characterization of new microsatellite 
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markers from Pacific abalone in Korea and they reported microsatellite markers could detect significant 
differentiation between one wild population and one hatchery population [17,18]. 

In the present study, we used six microsatellite markers to analyze the genetic diversity and 
relationships within and between cultured populations of H. discus hannai from different regions in 
Korea. This study will provide useful data for the effective monitoring and management of abalone 
populations as well as for the implementation of a stock enhancement program. 

2. Materials and Methods 

2.1. Sample Collection and DNA Extraction 

For the analysis, 223 Pacific abalones (H. discus hannai) were collected from five coastal locations in 
Korea in 2003 and 2004 (Figure 1). To cover the main coastal areas in Korea, five hatchery populations 
were sampled from three different areas (Table 1). The samples were obtained from hatchery-reared 
populations used as broodstock for artificial reproduction in Ulsan, in the eastern coastal areas (US; 5 1 
individuals); Namhae (NH; 49 individuals) and Wando (WD; 41 individuals) in the southern coastal 
areas; and Taean (TA; 39 individuals) in the western coastal areas. The samples for one hatchery-reared 
population were acquired from a population of hatchery-reared offspring in Geoje (GJ; 43 individuals) in 
the southern coastal areas. The WD hatchery managed a large broodstock composed of wild-caught and 
hatchery-produced adult Pacific abalones, the GJ hatchery kept no broodstock but has used 
hatchery-produced adults as parents for reproduction, and the others had small broodstocks composed of 
hatchery-produced adults. Although the original ancestors of the hatchery abalone were collected locally 
and the hatchery populations had been reared continuously, no detailed records of their founding and 
maintenance were available. 

Figure 1. Sampling sites and abbreviated names of five hatchery-reared Pacific abalone 
(H. discus hannai) samples (A) in Korea. The abbreviations are as follow: US (Ulsan), 
GJ (Geoje), NH (Namhae), WD (Wando) and TA (Taean). 
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Table 1. Collection details for five hatchery-reared populations of Haliotis discus hannai. 
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Sample Size 



Collection 



Date 



Ulsan population (US) Eastern Area (Ulsan); 35° 32' N, 129° 25' E 

Geoje population (GJ) Southern Area (Geoje); 34° 53' N, 128° 41' E 

Namhae population (NH) Southern Area (Namhae); 34° 42' N, 128° 0' E 

Wando population (WD) Southern Area (Wando); 34° 23 ' N, 126° 5 1 ' E 

Taean population (TA) Western Area (Taean); 36° 40' N, 126° 16' E 



41 



49 



43 



39 



51 



February. 2004 
June. 2003 



May. 2003 
March. 2004 



June. 2003 



Mantle musculature clips were preserved in 99.9% ethanol before being transported to the lab. Total 
DNA was extracted using the automated DNA extraction system MagExtractor MFX-2100 (Toyobo) 
with a MagExtractor-Genomic DNA Purification Kit (Toyobo, Osaka, Japan). Extracted genomic DNA 
was stored at -20 °C prior to polymerase chain reaction (PCR) analysis. 

2.2. Microsatellite Genotyping 

In total, 223 Korean Pacific abalones from the five populations were genotyped using six Pacific 
abalone microsatellite loci. Hdhl321, Hdh513, Hdh512 and Hdhl45 developed for Haliotis discus 
hannai [15] and Hddll4B and Hdd229 developed for Haliotis discus discus [14] were used to amplify 
alleles by PCR. The 5 '-end of the forward primer of each set of primers was labeled with fluorescent dye 
(6-FAM, HEX, or NED; Applied Biosystems, Foster City, CA, USA). PCR amplification of the six 
microsatellite loci was performed in 10-jxL volumes containing 0.25 U Taq DNA polymerase, 
10x ExTaq buffer, 2 mM dNTP mixture (Takara, Shiga, Japan), 2 jjM of each primer set and 
approximately 10 to 50 ng template DNA using a PTC-0220 DNA Engine Dyad Peltier thermal cycler 
(MJ Research, Inc., Waltham, MA, USA). PCR conditions included an initial denaturation at 95 °C for 
11 min, followed by 35 cycles of denaturation at 94 °C for 1 min, annealing for 1 min at each primer 
temperature listed in Table 1, and extension at 72 °C for 1 min, with a final extension at 72 °C for 5 min. 

For genotyping, 1 (iL of PCR product was added to 9 jxL of a reaction containing formamide (Hi-Di 
Formamide, Applied Biosystems, Warrington, UK) and the GeneScan 400HD [ROX] size standard 
(ABI PRISM, Applied Biosystems, CA, USA), denatured at 95 °C for 2 min, and immediately chilled on 
ice. Fragment analysis of the reaction products was performed using an ABI 3130 Genetic Analyzer 
(Applied Biosystems) and GeneMapper software (ver. 4.0; Applied Biosystems). To improve accuracy 
when determining allele sizes, a control DNA sample was included in each set of samples for each run. 

2.3. Data Analysis 

Statistical genetic analyses were conducted for five populations of H. discus hannai. The data were 
tested for PCR errors due to null alleles, stuttering, and allele dropout using MICRO-CHECKER 
ver. 2.2.3 (1000 randomizations) [21]. Scoring and human error were estimated by duplicate analyses. 
The data were also tested whether one or more microsatellites were under selection using LOSITAN 
software [22,23], which is a selection detection workbench based on a well evaluated F S T-outlier 
detection method. Two genetic diversity parameters, the number of alleles per locus (Na) and the number 
of unique alleles (U), were determined for each local sample at each locus, using the program Genepop 
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ver. 4.0 (http://kimura.univ-montp2.fr/~rousset/Genepop.htm). Allelic richness (/4r) was corrected for 
the smallest sample size (n = 39), using the rarefaction method of FSTAT ver. 2.9.3.2 [24]. Using FSTAT, 
allelic richness can be directly compared among populations, regardless of sample 
size [25]. For the analysis of molecular variance (AMOVA) [26], components of variance within and 
between populations based on the infinite allele model (IAM) were estimated using the program 
Arlequin ver. 3.0 [27]. The significance of AMOVA components was tested using 1,000 permutations. 
To estimate genetic heterozygosity among the entire set of pairwise population samples, unbiased 
expected (H e ) and observed (H 0 ) heterozygosity values [28] were calculated. Overall inbreeding 
coefficients (F IS ) [29] for each population and locus were also estimated; F IS can measure deviations 
from the Hardy- Weinberg equilibrium (HWE) within a population. Deviations from the HWE were 
tested using probability tests or exact tests by the Markov-chain procedure of Arlequin. 

The extent of population subdivision was examined by calculating global multilocus Fst values [29] 
and Rst values (1000 permutations) [30]. The index of pairwise F ST based on an infinite alleles model 
(IAM) was estimated and Rst values based on a stepwise mutation model (SMM) were calculated using 
Arlequin. The index i?sT incorporates the correlation of the weighted mean allele size, expressed as the 
number of tandem repeats. Significance levels were adjusted for multiple tests using the sequential 
Bonferroni correction technique [31]. 

The genetic distance between populations was estimated on the basis of the chord distance, Dce [32]. 
Dqe is one of the most efficient distance measures in obtaining correct tree topology from allele 
frequency data [33]. To examine genetic relationships between populations, a phylogenetic tree was 
constructed on the basis of pairwise genetic distances for all samples using the UPGMA method of 
the program POPULATIONS ver. 1.2.30 (http://bioinformatics.org/~tryphon/populations/). Bootstrap 
values were calculated using 1,000 replicates. The UPGMA tree was visualized using the Tree Explorer 
program (http://evolgen.biol.metro-u.ac.jp.TE/). 

Relationships between geographical populations were assessed using principal component analysis 
(PCA) of GenAlEx 6.3 (http://www.anu.edu.au/BoZo/GenAlEx/), which was based on the covariance 
matrix of gene frequencies. 

Because hatchery populations are often subjected to founder effects and bottlenecks that result in 
lower genetic diversity, Bottleneck software ver. 1.2.02 [34] was employed to test the bottleneck 
hypothesis under a two-phased model of mutation (TPM). This method can be used for testing the 
departure from mutation-drift equilibrium based on heterozygosity excess or deficiency. 

3. Results 

3.1. Genetic Variability 

Micro-Checker analysis did not detect any allele scoring errors caused by stuttering or large allele 
dropout, but it indicated that five of six loci could carry null alleles (Hddll4B was not affected; 
P < 0.05). However, all six loci were used in this study because no null alleles affected all populations. 

Lositan analysis did not detect any outlier loci that had excessively high or low Fst compared to neutral 
expectations, indicating all microsatellites used were not candidates for being subject to selection. 
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Estimates of genetic variability for the five Pacific abalone populations are summarized in Table 2. 
These results suggest that all of the microsatellite loci were polymorphic, with some level of 
polymorphism, including large differences in the number of alleles, in all of the studied populations. Na 
and A R ranged from 20.17 in the TA population to 26.67 in NH and from 20.17 in TA to 24.78 in NH, 
respectively, with the degree of variability differing among the six loci. The mean expected and observed 
heterozygosities per locus ranged from 0.885 in WD to 0.919 in NH and from 0.752 in NH to 0.838 in 
WD, respectively. However, despite these differences in genetic diversity, no significant decrease in 
genetic variability was found among the population samples (Kruskal-Wallis test, P > 0.05). In total, 56 
alleles were found to be unique to a single population. The greatest number of unique alleles was 
detected in NH (20), while TA had only seven unique alleles. The inbreeding coefficient Fis of the 
hatchery populations ranged from 0.054 to 0.184 (meanFis = 0.121 ± 0.056). 

The observed genotype frequencies were tested for agreement with HWE (Table 2). Among the 30 
population-locus cases, in 18 cases the observed genotype distribution was generally in accordance with 
Hardy- Weinberg proportions; however, 12 (40%) departed significantly from HWE, even after 
adjustment of P values using the sequential Bonferroni method (P < 0.01), and these disequilibriums 
(DHWE) were detected at all loci except Hdhll4B. Ten DHWE cases resulted from a deficiency of 
heterozygotes and two from an excess of heterozygotes (Table 2). Eleven of the 30 Fis values at the six 
loci estimated for the five populations were significantly different from zero (P < 0.01). Significant 
deviations were not evenly distributed among samples or loci, nor were they associated with a particular 
locus or sample. DHWE was observed at one locus in WD, but at four loci in GJ. 

In total, 233 different alleles were observed across all loci and all samples, ranging from 15 for the 
Hdhl54 locus to 64 for the Hdh513 locus. No population had a diagnostic allele. In the five populations, 
64.0 to 76.6% of the alleles were rare alleles with a frequency <5%. They were detected at most loci and 
were not associated with a particular locus in either population. 

In testing the departure from mutation-drift equilibrium based on heterozygosity excess or deficiency 
for five populations, bottleneck analysis was conducted using Bottleneck software under the TPM of 
microsatellites. No population displayed significant heterozygosity excess (P > 0.05) through the sign 
test, standardized differences test and Wilcoxon sign rank test, suggesting that five hatchery populations 
have not experienced a recent bottleneck. 

Table 2. Allelic variability observed at six microsatellite loci in five hatchery-reared 
populations of Haliotis discus hannai. 



Population (No.) 



Microsatellite Loci 

Hdhl321 Hddll4B Hdd229 Hdh513 Hdh512 Hdhl45 Mean 



Ulsan (51) 



Na 
A R 
S 

u 



29 26 22 39 31 10 26.17 

27.59 24.07 20.40 34.45 28.70 9.23 24.07 
264-380 152-256 168-218 422.000 70-200 130-154 




2 2 0 2 4 1 1.83 

0.961 0.943 0.913 0.967 0.964 0.750 0.916 
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Table 2. Cont. 



Population (No.) 




Microsatellite Loci 
Hdhl321 Hddll4B 


TTHrl990 




1 1 1 1 1 1 . ~ i _ 




iviean 




TT 


1.000 


0.902 


n. £/n 


n. 89/i 


1 f\(\(\ 

1 .uuu 


n ^/io 
u. j^y 


f» 89H 
U.oZU 




J7„- 
r IS 


-0.041 


0.044 


A OAT* 


n. 1 aq* 


U.Uj / 


u.z /u 


U.1UO 




r> 
r 


0.934 


0.727 


u.uuu 


U.Ujo 


n 9 1 1 

U.Zl 1 


u.uuu 




Geoje (43) 


N A 


9Q 


9^ 

Z-J 


21 


31 


23 


1 A 

10 


22.83 




Ar 


ZO.UZ 


99 

ZZ.U J 


20.42 


29.83 


22.52 


9.71 


22.09 




s 

kJ. 


Z04-JOO 


1 ^9 9^9 
1 jZ-ZjZ 


172-226 


190-422 


88-148 


114-154 






u 


A 

U 


1 


2 


3 


2 


1 


1.50 




H. 


Kl.yjy 


fi 8Q8 


0.905 


0.958 


0.950 


0.693 


0.894 




H„ 


n. s 1 a 

U.014 




0.674 


0.744 


0.953 


0.488 


0.767 




F rc 
1 lb 


U. 1 j j 


-V.VJ 1 


0.257* 


0.225 


-0.003 


0.298* 


0.143 




P 


u.uuu 




0.000* 


0.011 


0.001* 


0.000* 




Namhae (49) 


N A 


32 


29 


O 1 
Z 1 


/in 
4U 


z / 


1 1 
1 1 


26.67 




Ar 


98 Q9 


97 fi4 
Z / .UH 


20.16 


36.38 


25.36 


1 A O ^> 

10.83 


24.78 




s 


988 1Q8 
Zoo— jyo 


1 'sA 9^8 
1 j4— Zjo 


164-218 


174-422 


74-142 


106-152 

J. \J\J l ^/ Z- 






u 


f. 


1 

1 


1 


8 


1 


3 


3.33 




H. 


0.925 


0.951 


0.945 


0.968 


0.950 


0.777 


0.919 




H 0 


0.571 


0.918 


0.857 


0.735 


0.878 


0.551 


0.752 




Fis 


0.385* 


0.034 


0.094 


0.243* 


0.077 


0.293* 


0.184 




P 


0.000* 


0.485 


0.128 


0.001* 


0.095 


0.000* 




Wando (41) 


Na 


26 


24 


18 


36 


20 


5 


21.50 




Ar 


25.60 


23.46 


17.56 


35.16 


19.89 


5.00 


21.11 




S 


988-^78 


1 ^9-970 


168-224 


182-422 


74-132 


130-140 






u 


2 


4 


1 


2 


0 


0 


1.50 






A A A C 

0.945 


A A'? A 

0.930 


ft 87Q 

u.o / y 




ft Q^4 
\j .y 




0.885 




T T 


0.951 


0.927 


0.683 


0.976 


A A f 1 

0.951 


0.537 


0.838 




Fis 


-0.007 


0.004 


0.225* 


-0.007 


-0.019 


0.178 


0.054 




P 


0.157 


0.764 


0.000* 


0.178 


0.275 


0.302 




Taean (39) 


Na 


28 


19 


18 


31 


19 


6 


20.17 




Ar 


28.00 


19.00 


18.00 


31.00 


19.00 


6.00 


20.17 




S 


266-384 


152-260 


172-218 


188-422 


90-134 


128-140 






u 


3 


2 


0 


2 


0 


0 


1.17 




H e 


0.945 


0.907 


0.925 


0.964 


0.922 


0.781 


0.907 




H 0 


0.974 


0.897 


0.462 


0.897 


0.974 


0.641 


0.808 




Fis 


-0.031 


0.010 


0.505* 


0.070 


-0.058* 


0.181 


0.111 




P 


0.996 


0.581 


0.000* 


0.170 


0.000* 


0.025 
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Table 2. Cont. 



Population (No.) 



Microsatellite Loci 

Hdhl321 Hddll4B Hdd229 Hdh513 Hdh512 Hdhl45 Mean 



Mean all populations 




N A 
A R 
U 
H e 
Ho 



28.80 

27.63 

2.60 

0.951 

0.919 



24.20 

23.12 

2.00 

0.926 

0.915 



20.00 

19.31 

0.80 

0.909 

0.607 



35.40 

33.36 

3.40 

0.965 

0.835 



24.00 

23.10 

1.40 

0.944 

0.951 



8.40 

8.15 

1.00 

0.730 

0.553 



23.47 

22.44 

1.87 

0.904 

0.797 



Number of samples (No.), number of alleles per locus (N A ), allellic richness (^4 R ), size in bp of alleles (S), 
number of unique alleles (U), expected heterozygosity (He), observed heterozygosity (Ho), inbreeding 
coefficient (F IS ), and probability of significant deviation from Hardy- Weinberg equilibrium (P) are given 
for each population and locus. Wide significance levels were applied using the sequential Bonferroni technique 
(k= 6) [31]. * Significant at P < 0.01. Calculations assume that individuals with one microsatellite band are 
homozygous for the allele. Number in parentheses below F ls indicates the probability of significant 
heterozygosity excess or deficit. 

3.2. Genetic Differentiation among Populations 

Low but significant genetic differentiation (overall F S t = 0.009, P<0.01) was observed among 
populations, indicating that the level of genetic heterogeneity among them was low. Table 3 shows the 
pattern of genetic differentiation among populations observed by comparing Dqe and mean pairwise Fst 
and Rst values using composite allele frequency data. Values of pairwise Fst among populations were 
significantly different from zero in all pairwise comparisons except one (all P < 0.01 after sequential 
Bonferroni correction, except for P = 0.086 between NH and US). On the other hand, the pairwise Rst 
among populations was significantly different from zero in only five pairwise comparisons. Estimates of 
F ST (range 0.0040 to 0.0176) were generally lower than R ST values (range -0.0082 to 0.1016). The D CE 
values were similar to the F$ T and Rst values, in the range of 0.018 to 0.026. The AMOVA from the 
distance matrices for 233 individuals permitted partitioning of the overall variation. The amount of 
variation attributable to differences between individuals within populations and within-individual 
differences were large (12.09 and 87.06, respectively; P < 0.05), while only a low but significant 
proportion of variation (0.85%; P < 0.05) was detected among populations. The UPGMA tree and the 
PCA scatter plot constructed on the basis of the Dqe values showed that no population formed a separate 
cluster (Figures 2 and 3). The US population from the eastern coast and the NH population from the 
southern coast were the most closely related populations, while GJ from the southern coast was the most 
isolated from the others, the same groupings found with the PCA scatter plot analysis. Together, pairwise 
Fst and Rst tests, the UPGMA tree, and PCA on the five cultured populations of H. discus hannai did not 
show a consistent relationship between geographic and genetic distances. 
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Table 3. Dce distance (below the diagonal) and mean F ST estimates and Rst estimates 
(above the diagonal) between each pair of five hatchery-reared populations of Haliotis 
discus hannai. 

Population Ulsan Geoje Namhae Wando Taean 



Ulsan 




0.0068* 
(0.0744*) 


0.0050* 
(0.0341 NS ) 


0.0088* (0.0198 NS ) 


0.0091* (0.0053 NS ) 


Geoje 


0.019 




0.01 15 NS 
(0.0289*) 


0.0124* (0.023 1 NS ) 


0.0176* (0.1016*) 


Namhae 


0.018 


0.024 




0.0068* (-0.0082 NS ) 


0.0098* (0.0524*) 


Wando 


0.019 


0.021 


0.020 




0.0172* (0.0332*) 


Taean 


0.018 


0.026 


0.023 


0.020 





D CE distance [32] and pairwise F ST and R ST [29,30] are measures of genetic distance and genetic differentiation 
between populations, respectively. Number in parenthesis indicates R$j. Wide significance levels were applied 
using the sequential Bonferroni technique (k = 10) [31]. * Significant at P < 0.01. NS is nonsignificant after 
sequential Bonferroni correction. 

Figure 2. UPGMA dendrogram showing the phylogenic relationships among five cultured 
populations of the Pacific abalone. The abbreviations are as follow: US (Ulsan), GJ (Geoje), 
NH (Namhae), WD (Wando) and TA (Taean). 

GJ 

49| US 

40[1 NH 

I TA 

48 

I WD 

0.05 



Figure 3. Principal components analysis, plotting the relationships of the studied five 
H. discus hannai cultured populations based on Dqe genetic distances [21]. The abbreviations 
are as follow: US (Ulsan), GJ (Geoje), NH (Namhae), WD (Wando) and TA (Taean). 
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U 

Oh 
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4. Discussion 

Although the two loci, Hddll4B and Hdd229, used in this study were originally described and 
designed for a congener (H. discus discus) [14], they all worked satisfactorily for//, discus hannai. This 
remarkable level of cross-species microsatellite conservation is a clear indication of the close molecular 
relatedness of the two species. Similar results were observed in the microsatellite analysis of H discus 
discus [17,35]. 

The genetic variability of the six microsatellite loci proved to be extensive. Genetic variability in the 
five hatchery populations (average Na = 23.47, average H e = 0.904) was significantly higher than 
previously measurements from Japanese and Chinese hatcheries (average Na = 5.43 to 8.7, average 
H e = 0.575 to 0.774) [6,20]. However, Li et al. [6] reported that for the Japanese natural Pacific abalone 
populations, the average number of alleles per locus and heterozygosity was 22.4 and 0.793, 
respectively. Similar genetic variability were also detected for eight Korean natural Pacific abalone 
populations (average Na = 28.13, average He = 0.916; An et ah, data not published). Most of 
microsatellite markers used in above studies were same; then, can be compared on an equal basis. These 
values are similar to our results, indicating that a high level of genetic variation has been preserved 
during the domestication of these five cultured Pacific abalone populations in Korea. The maintenance 
of genetic diversity in hatchery populations has been reported for other marine organisms. For example, 
no differences in genetic diversity could be detected between wild and hatchery populations of Pacific 
oysters in Australia and China [36,37]. The use of hundreds of parental males and females from different 
farms for spawning might contribute to the maintenance of genetic diversity in these animals. 

Significant deviations from HWE were observed in 12 of the 30 population-locus cases. Ten cases 
had lower heterozygosity than predicted by HWE. Heterozygote deficiency has been reported for many 
other marine invertebrates [38,39]. However, in many cases, its cause remains unknown [40]. In the 
hatchery populations, heterozygote deficiency is commonly caused by a limited number of founders or 
founder effects. The observed Fis in this study can attribute to the founder effect. The limited dispersal 
nature of abalones together with the sweepstake style of reproductive success may easily lead to such a 
non-random mating of population [35]. The heterozygote deficiency observed in hatchery populations is 
probably caused by unequal contributions of a limited number of pair mating to the next generation. 
Alternatively, the null alleles observed at the microsatellite loci present a likely explanation for the 
deviation [41]. In our study, null alleles were observed at five of the six loci, and no significant deviation 
from HWE was observed at the Hddll4B locus, which did not have any null alleles. In contrast, two of 
the 12 significant HWE deviations were caused by an excess of heterozygotes. Recent population 
bottlenecks, when a population of interest experiences a reduction in size, sometimes also reduce the 
number of homozygotes [42,43]. However, bottleneck analysis showed that five hatchery populations 
have not experienced a recent bottleneck; thus, the possibility of recent bottleneck seems unlikely. 
Alternatively, it is possible to ascribe the heterozygote excess in cultured stocks to an overdominance 
phenomenon, which reduces the survival of homozygotes [44]. Li et al. [6] and Hara and Sekino [19] 
observed, on the basis of microsatellite analysis, heterozygote excesses in 8 of 10 and 13 of 16 
significant HWE deviations in Pacific abalone hatchery stocks, respectively. 

Genetic differentiation detected by Fst and Rsi values were low but significant beteen most pairs of 
hatchery populations. Although we cannot trace the origins of the significant differences between the 
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hatchery populations with no records from the farms, this differentiation observed among cultured 
populations is likely be due to a low number of founding individuals in hatcheries, which increases the 
effect of genetic drift, as has previously been reported for other hatchery-reared species including 
abalone [6,19,45^17]. In addition, different selection procedures in hatchery practice may have also led 
to changes in the genetic composition of hatchery populations [20]. The AMOVA of all six 
microsatellites revealed that 0.85% (Fst = 0.009, P = 0.040) of the genetic variance occurred among the 
hatchery samples, with the remainder of the variance occurring between individuals within populations 
and within individuals. To our knowledge, this is the lowest among-population variance component 
reported for this species. Previously reported F S t values for hatchery Pacific abalone (F ST = 0.243 to 

0. 427 at six polymorphic loci for three hatchery populations in Japan; Fst = 0.056 to 0.259 at nine 
polymorphic loci for nine hatchery populations in Japan; Fst = 0.0546 to 0.0976 at seven polymorphic 
loci for five hatchery populations in China) [6,19,20] all exceeded our estimated value. The UPGMA 
tree and the PCA scatter plot constructed using Dqe genetic distances did not indicate a regional structure, 

1. e., individuals from nearby regions were not grouped together (Figures 2 and 3). Although differences 
in the number of polymorphic loci and the number of populations examined in a study will influence F$t 
values and some element of hatchery selection may have also led to changes in the genetic composition 
of the cultured populations, the maintenance of genetic diversity, the relatively low Fst value and the 
lack of a relationship between the geographic and genetic distances suggest that some gene flow has 
occurred among populations during the extensive aquaculture of abalone, though the possibility of 
widespread exchanges of stocks and eggs between hatcheries by local farmers cannot be excluded. 

5. Conclusions 

In conclusion, an understanding of the genetic structure and diversity among marine resources, 
especially for species for which artificial stocks are produced by hatchery-produced seed for natural 
resource restoration, are critical for the establishment of suitable guidelines for resource management 
and selective breeding. No detailed information is available to date on the genetic diversity of hatchery 
populations of Pacific abalone in Korea. In this study, we reported that relatively high genetic variability 
and significant but minor genetic differentiation were detected among the hatchery populations of 
Pacific abalone in Korea based on analysis of six microsatellite loci. The lack of a relationship between 
the geographic and genetic distances implied extensive abalone aquaculture. The allelic composition and 
diversity of Pacific abalone should be carefully considered regarding the seedling and stocking practice 
of hatcheries in order to conserve the genetic diversity of the natural population. Such information will 
be of assistance in the genetic management of fisheries and the successful implementation of 
stock-enhancement programs. 
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